Submitted to the Astrophysical Journal 

Preprint typeset using I^TgX style emulateapj v. 05/04/06 



oo 
O 
O 
(N 

-)— > 
o 

O 

(N 
(N 

o 



(N 
> 

rn 

O 
oo 
O 



X 
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ABSTRACT 

We present a two-dimensional analysis of the bright nearby galaxy group NGC 5044 using the currently 
available Chandra and XMM data. In the inner 1 kpc a pair of cavities are evident together with a set of bright 
X-ray filaments. If the cavities are interpreted as gas displaced by relativistic plasma inflated by an AGN, even 
in the absence of extended 1 .4 GHz emission, this would be consistent with a recent outburst as also indicated 
by the extent of dust and Ha emission. The soft X-ray filaments coincident with Ha and dust emission are 
cooler than the ones which do not correlate with optical and infrared emission. We suggest that dust-aided 
cooling contributes to form warm (T ^ 10"* K) gas, emitting Ha radiation. At 31 kpc and 67 kpc a pair of cold 
fronts are present, indicative of sloshing due to a dynamical perturbation caused by accretion of a less massive 
group, also suggested by the peculiar velocity of the brightest galaxy NGC 5044 with respect to the mean group 
velocity. 

Subject headings: cooling flows — galaxies: clusters: general — galaxies: clusters: individual 
(NGC 5044) — X-rays: galaxies: clusters 



1. INTRODUCTION 

The current X-ray observatories, Chandra and XMM, have 
revolutionized our understanding of the cores of relaxed 
galaxies, groups and clusters, which show a highly peaked X- 
ray emission from a hot interstellar medium whose radiative 
cooling time is less than 1 Gyr (for recent determinations of 
cooling times in these objects see for example Voigt & Fabian 
2004; Sanderson et al. 2006; Jetha et al. 2007). In the ab- 
sence of heating, a cooling flow is established, in which the 
gas cools, condenses and flows toward the center, accreting 
onto the central galaxy (Fabian 1994). However the mass 
sink for all this supposedly cooling and condensing gas has 
never been found (Donahue & Voit 2004). X-ray observa- 
tions with Chandra and XMM have established that there is 
little evidence for emission from gas cooling below ~ T^u-/3: 
just when gas should be cooling most rapidly it appears not to 
be cooling at all (see the recent review by Peterson & Fabian 
2006, and references therein). 

A compensating heat source must therefore resupply the 
radiative losses, and many possibilities have been proposed, 
including thermal conduction (e.g., Narayan & Medvedev 
2001), energy released by mergers (e.g., Motl et al. 2004) 
or by supernovae (e.g.. Silk et al. 1986). However, feedback 
from the central AGN has become the most appealing solution 
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to the problem (see the recent review by McNamara & Nulsen 
2007, and references therein). There is, in fact, clear obser- 
vational evidence for AGN heating as the majority of bright- 
est central galaxies of clusters and groups host a radio loud 
AGN (e.g.. Burns 1990; Best et al. 2007) and, following the 
launch of Chandra, in an increasing number of objects such 
disturbances as shocks, ripples and cavities have been found 
in the central atmospheres of clusters, groups and elliptical 
galaxies (e.g., Fabian et al. 2006; Birzan et al. 2004; Dunn & 
Fabian 2006; Forman et al. 2005; Vrtilek et al. 2002; Jetha 
et al. 2007; Croston et al. 2008; Finoguenov & Jones 2001; 
Allen et al. 2006). The cavities, which appear as X-ray sur- 
face brightness depressions, have been interpreted as bubbles 
of low density relativistic plasma inflated by radio jets, dis- 
placing the thermal gas causing PdV heating (e.g., Churazov 
et al. 2002), although Mathews & Brighenti (2008) pointed 
out that cavity formation contributes thermal energy that may 
offset radiative cooling only by injecting ultra-hot but non rel- 
ativistic gas; X-ray cavities formed solely from relativistic gas 
have a global cooling effect. In most cases, the energy intro- 
duced by the AGN is more than sufficient to counteract puta- 
tive cooling flows (McNamara & Nulsen 2007), although the 
physical process of the coupling of the feedback energy with 
the ambient medium is not well understood. While many of 
the observed X-ray cavities are filled with plasma emitting in 
the radio at 1 .4 GHz, some are undetected at this frequency 
and have been referred to as "ghost cavities". These may re- 
sult from the aging of the relativistic particle population and 
be the signature of a previous AGN outburst. Observations at 
low radio frequency of several rich clusters with ghost cavities 
show that they are indeed filled with relativistic plasma (e.g., 
A 2597, Clarke et al. 2005), although even low frequency 
emission cannot be easily detected. 

However, some net cooling is probably occurring, as the 
galaxies at the center of cool cores show properties not shared 
by typical elliptical galaxies (Crawford 2004): 1) the pres- 
ence of strong, low-ionization emission line nebulae both in 
clusters (e.g., Crawford et al. 1999; Conselice et al. 2001) 
and in groups (e.g. Macchetto et al. 1996); 2) the galaxies 
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at the center of clusters that display these nebular emission 
show a component of excess ultraviolet/blue continuum asso- 
ciated with young stars (e.g., Cardiel et al. 1998; Crawford 
et al. 1999). This excess blue light is again extended on scales 
of several kpc and there is evidence that some (spatially ex- 
tended) star formation has been triggered by interaction with 
the radio source (McNamara 2004) 3) CO emission lines have 
been detected in several cool core clusters at millimeter wave- 
lengths (e.g., Edge 2001; Salome & Combes 2003), implying 
the presence of a substantial amount of warm molecular gas 
(10^""^ M0) within a 50 kpc radius of the central galaxy. 
Less massive warmer molecular regions have been observed 
in H2 lines, which are often spatially associated with the stel- 
lar UV and Ha emission extending over ^ 20 kpc (Jaffe et al. 
2005). 

The origin and excitation mechanism of the nebular Ha 
emission have been much debated and they are still poorly 
known. The nebulae require a constant and distributed heating 
source (e.g. Johnstone & Fabian 1988) which could be stellar 
(e.g. massive OB stars, Allen et al. 1992) or the ICM (e.g. 
conduction. Sparks et al. 1989). Deep Chandra X-ray imag- 
ing of Perseus (Fabian et al. 2003b) and M87 (Sparks et al. 
2004) have renewed interest in conduction because of the spa- 
tial coincidence of filamentary Ha and soft X-ray emission, 
which can be due to conduction and mixing of the cold gas 
with the ICM. The disposition of some of the optical fila- 
ments in Perseus strongly suggests that the filaments are due 
to the buoyant radio bubbles drawing out the cold gas from a 
central reservoir (Fabian et al. 2003b) and it reveals the pos- 
sibility of tracing the flow of the rising bubbles (Hatch et al. 
2006). But the question is still open if galaxy-galaxy interac- 
tions are required to stimulate the central optical nebulae and 
if the optical filaments are trails of galaxies punching through 
a molecular hydrogen reservoir as suggested by the examples 
in Wilman et al. (2006). The excitement caused by the shift 
in the cooling flow paradigm is stimulating also a renewed 
theoretical interest in the quest for the origin and excitation 
of molecular (Ferland et al. 2008) and Ha emission (Nipoti 
& Binney 2004; Pope et al. 2008; Revaz et al. 2008). In par- 
ticular Nipoti & Binney (2004) have shown how the Ha neb- 
ulae can be stable against thermal evaporation only for the 
conditions of temperature and pressure found in cool cores, 
explaining their association with that type of environment. 

High resolution Chandra observations have also revealed 
another interesting and unexpected feature in the ICM of re- 
laxed clusters: the presence of cold fronts. Many clusters have 
been found to exhibit sharp arc-shaped jumps in their gas den- 
sity and temperature which, unlike shock fronts, have the gas 
on the dense side cooler, so the pressure is continuous across 
the front (e.g., Markevitch et al. 2000; VikhHnin et al. 2001). 
In merging systems they were interpreted as contact disconti- 
nuities between gases from different sub-clusters (Markevitch 
et al. 2000). However cold fronts are present in the centers 
of many, if not most, relaxed clusters with cool cores (e.g., 
Mazzotta et al. 2001; Markevitch et al. 2003; Dupke & White 
2003; Ascasibar & Markevitch 2006; Ghizzardi et al. 2007). 
The widely accepted scenario is that cold fronts in these sys- 
tems are due to sloshing of the cool gas in the central gravita- 
tional potential, which is set off by minor mergers/accretions; 
the only necessary condition is a steep entropy profile as ob- 
served in relaxed clusters (Ascasibar & Markevitch 2006). 

Whereas a growing number of clusters and elliptical galax- 
ies have deep enough multi-wavelength data (X-rays, radio 
and optical) to study the rich phenomenology of cool cores 



in a detailed spatially resolved fashion, only an handful of 
groups with such coverage exists, as for example the GEMS 
objects (Osmond & Ponman 2004) presented in Rasmussen 
& Ponman (2007), and therefore, "unfortunately, AGN heat- 
ing is not as well studied in groups as in clusters" (McNa- 
mara & Nulsen 2007). Examination of AGN feedback at the 
mass scale of groups is valuable because, although the scale 
of outbursts in groups is less energetic and often on a smaller 
spatial scale than in clusters, the impact can be even more 
dramatic than in rich clusters due to the shallower group po- 
tential. Statistical studies examining the impact of AGN on 
groups are starting to address the points raised above (Cros- 
ton et al. 2005; Jetha et al. 2007). Cold fronts have not been 
investigated in detail so far in relaxed groups of galaxies. 

In this paper we show with more detail the currently avail- 
able Chandra and XMM data for the galaxy group NGC 5044, 
one of the brightest groups in X-rays: it is one of the only 5 
objects with kr < 2 keV in the HIFLUGCS sample (Reiprich 
& Bohringer 2002). Our estimate, using the new Chandra 
and XMM data, for the bolometric (0.1-100 keV) X-ray lu- 
minosity within rsoo = 443hjQ kpc (Gastaldello et al. 2007b) 
is 1.05 ±0.06 X lQ*^hjQ erg s~\ in good agreement with the 
ROSAT estimate presented in Reiprich & Bohringer (2002) 
of 2.46 X lO^^h'^Q within 560h'^Q kpc. In our previous analy- 
sis, mainly focused on radial, azimuthally averaged properties 
of the system (Buote et al. 2003a,b), interesting features like 
holes and filaments within 10 kpc, a sharp edge, resembling 
a cold front, at 67 kpc, and unusually low iron abundances at 
large radii (Buote et al. 2004) were already pointed out. An 
Ha nebula is present in the core of NGC 5044 showing an 
extended filamentary structure (Goudfrooij et al. 1994a; Caon 
et al. 2000). NGC 5044 is also remarkable because Spitzer 
data show extended cold dust emitting at 70 /im (Temi et al. 
2007a). Moreover, Temi et al. (2007b) showed extended 8 ^m 
excess (likely arising from PAH, polycyclic aromatic hydro- 
carbon, molecules) extending out to several kpc and spatially 
coincident with the Ha emitting nebulosity and the brightest 
soft X-ray emission. As proposed in Temi et al. (2007b), cur- 
rent evidence is consistent with an internal origin of this dust, 
which has been buoyantly transported from the galactic core 
out to several kpc into the hot X-ray emitting gas following 
an AGN outburst. 

All distance-dependent quantities have been computed as- 
suming //()= 70 km s"' Mpc"', Qm= 0.3 and fiA= 0.7. At the 
redshift of z = 0.009 1' corresponds to 11.1 kpc. All the errors 
quoted are at the 68% confidence limit. 

2. OBSERVATIONS AND DATA PREPARATION 

NGC 5044 has been observed by Chandra with the ACIS-S 
configuration on September 1999 (obsID 798) for 20.7 ks and 
by XMM on January 200 1 (obsID 0037950 1 1 ) with the EPIC 
MOS and pn cameras for 23 and 20 ks, respectively (see also 
Buote et al. 2003a,b). We focused mainly on the Chandra 
data because its PSF is much better suited for the high-spatial 
resolution study that constitutes most of the paper. We take 
advantage of the XMM data when we discuss the properties at 
larger radii (see section §3.2 and §4.2) where the resolution 
requirements are not so important and we can benefit from 
the larger accessible field. Here we provide just a brief de- 
scription of the data preparation, more details can be found in 
Gastaldello et al. (2007b). 

2.1. Chandra 
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Fig. 1. — Chandra 0.5-5 keV X-ray image of the inner 30x30 kpc of 
NGC 5044. The image has been processed to remove point sources, flat 
fielded with a 1.0 keV exposure map. The figure has been contour binned us- 
ing the algorithm of Sanders (2006) with a S/N of 10. Color bar units are total 
counts per pixel. Interesting features discussed in the text are highlighted. 
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Fig. 2. — Unsharp-masked Chandra image created by subtracting a 0.5-5 
keV intensity image smoothed by a ct = 10" Gaussian and one smoothed by 
a (T = 1" Gaussian. Areas in black show a deficit of counts, whereas areas in 
white show an excess of counts. Also shown are the angular sectors used for 
the surface brightness profiles shown in Fig. 3, labeled by name according to 
the scheme of Table 1 . 



The data were analyzed with the X-ray analysis packages 
CIAO 3.4 and Heasoft 6.4 in conjunction with the Chandra 
calibration database iCaldb) version 3.4.2. In order to en- 
sure the most up-to-date calibration, all data were reprocessed 
from the "level 1" events files, following the standard Chan- 
dra data-reduction threads'". We applied the standard cor- 
rections to take account of a time-dependent drift in the de- 
tector gain and charge transfer inefficiency, as implemented 
in the CIAO tools. From low surface brightness regions of 
the active chips we extracted a light-curve (5.0-10.0 keV) to 
identify and excise periods of enhanced background. The ob- 
servation was quiescent resulting in a final exposure time of 
20 ks. Point source detection was performed using the CIAO 
tool wavdetect and removed, so as not to contaminate the 
diffuse emission, using appropriate elliptical regions contain- 
ing 99% of their flux. We generated an image in the 0.5-5.0 
keV and a corresponding exposure map computed at an en- 
ergy of 1.0 keV. A zoom of the image in the inner region is 
presented in Fig.l and Fig. 10 and the the region covered by 
the ACIS-S3 chip is presented in Fig.4. 

2.2. XMM 

We generated cahbrated event files with SAS v7.1.0 using 
the tasks emchain and epchain. We considered only event pat- 
terns 0-12 for MOS and for pn, and the data were cleaned 
using the standard procedures for bright pixels and hot column 
removal and pn out-of-time correction. Periods of high back- 
ground due to soft protons were filtered as in Gastaldello et al. 
(2007b) resulting in a net exposure time of 22 ks for MOS 1 
and MOS2 and 17ks for pn. For each detector, we created im- 
ages in the 0.5-2 keV band with point sources detected using 
the task ewavelet and masked using circular regions of 25 
radius centered at the source position. We created exposure 
maps for each detector and we combined the MOS images 
into a single exposure-corrected image, smoothed on a scale 
of 10", shown in Fig. 5. 

'" http://cxc.harvard.edu/ciao/threads/index.html 



3. X-RAY IMAGES AND SURFACE BRIGHTNESS 
PROFILES 

3.1. The inner 10 kpc 

The presence of a disturbed morphology with filamentary 
structure is already evident from the raw Chandra image and 
further confirmed by using the contour binning technique 
of Sanders (2006) (See Fig.l). Two depressions in surface 
brightness and multiple filamentary structures, some of them 
connected to the presence of the cavities, have been high- 
lighted in Fig.l. To highlight the presence of structure we 
performed an "unsharp masking" of the Chandra image (e.g., 
Fabian et al. 2003a): the 0.5-5 keV exposure-map-corrected 
image was smoothed with Gaussians of width 1 and 10 arcsec 
and the two smoothed images were then subtracted, with the 
resulting image shown in Fig. 2. The southern cavity and the 
structure of filaments are clearly shown in the unsharp masked 
image. To quantify the observed structures we extracted from 
the exposure corrected image (the unsharp image was not used 
in the following analysis of the surface brightness profiles) an 
azimuthally averaged radial profile and profiles from selected 
angular sectors as depicted in Fig. 2 and listed in Table 1 (cen- 
tered on the X-ray surface brightness peak coincident with the 
optical center of the galaxy NGC 5044 as listed by NED and 
by Macchetto et al. 1996, RA 13:15:24.0 and DEC -16 23 09, 
J2000 coordinates). 

As it can be seen from Fig. 3, there is a complex network of 
filamentary structure of enhanced emission which highlights 
dramatically regions of comparatively lower surface bright- 
ness. Two filaments, brighter and extended (in particular com- 
pared to the azimuthally averaged profile, also shown in the 
right panel of Fig. 3), are evident in sector I (the SE filament) 
and sector III (the SW filament) with an embedded cavity in 
Sector II (also in this sector the emission is more extended 
compared to the average, i.e. there is an excess of counts for 
radii greater than 1'). The filament in sector III is brighter in 
the region 10"-30" but less extended, i.e. dimmer at r > 1', 
than the overall emission at PA 120-182. The emission in sec- 
tor IV can robustly be considered filamentary, with enhanced 
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Fig. 3. — Top panel: Surface brightness profiles in selected angular sectors listed in Table 1 and shown in Fig. 2. The azimuthally averaged surface brightness 
profile is plotted as a solid black line in both panels for ease of comparison. 
Bottom panel: Ratio of the surface brightness profiles of the above selected angular sectors over the azimuthally averaged one. 



TABLE I 



Sector 


PA 


Sector I 


120-141 


Sector II 


141-182 


Sector III 


182-212 


Sector IV 


212-244 


Sector V 


244-280 


Sector VI 


310-345 



Note. — Position an- 
gles (measured form the 
N direction) for ttie sur- 
face brightness profiles 
discussed in the text in 
§3.1. 



emission, even brighter than the other two filaments in the 
range 30"-50" but less extended, with an embedded depres- 
sion in surface brightness between 10"and 30", which stands 
out as a region of low emission, compared to the surround- 
ings, in Fig. 2. However this feature is at the same level of 
surface brightness of the azimuthally averaged profile. The 
cavity in sector VI is the real counter-part of the one in sector 
II. Although it stands out less clearly in the unsharp-mask im- 



age, it is clearly indicated in the binned images of Fig. 1 and in 
the striking similarity of the surface brightness profiles within 
50" as shown in Fig. 3. The northern cavity stands out less 
clearly in the unsharp masked image of Fig. 2 because of the 
dimmer regions surrounding it (the northern arms are not as 
bright as the southern ones) and because of the sharper drop 
in surface brightness at r > 50". There is also evidence of 
enhanced emission in sector V in the 20"-40" region, as high- 
lighted also in the unsharp masked image. 

3.2. The inner and outer surface brightness edges 

We now zoom out to investigate the outer regions in the 
Chandra image, as shown in Fig.4. A sharp edge in the sur- 
face brightness on the south-east (PA 120-160) is revealed 
in this image and further confirmed by a surface brightness 
profile in the same sector, as compared to a reference pro- 
file extracted in the PA 180-315 sector (right panel of Fig.4). 
Whereas the profile is smooth and with a continuous deriva- 
tive in the reference sector, the slope of the sector across the 
cold front is rapidly changing and the edge is at ~ 31 kpc 
(170"). 

Inspection of the combined MOS image shown in Fig. 5 
confirms the suggestion of Buote et al. (2003b) of the pres- 
ence of another surface brightness edge in the north-west di- 
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Fig. 4. — Left panel: Over-plotted on the Chandra image are the sector used to extract the surface brightness profile across the cold front and the reference 
sector used for comparison. 
Right panel: Surface brightness profile for the cold front and the reference sector A sharp drop is visible at ~ 170" form the center, i.e. 31 kpc. 
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Fig. 5. — Mosaic of the MOSl and MOS2 images smoothed on a 10" 
scale. The image has been divided by the summed exposure maps to correct 
for exposure variations. Point sources have been removed and filled using the 
CIAO task dmf ilth. Sectors used for the extraction of surface brightness 
profiles discussed in the text and presented in Fig. 6 are also shown. 

rection at ^ 6'. Plots of the surface brightness profiles across 
this sector and in other three sectors (one including the in- 
ner edge) for comparison are shown in Fig. 6. The inner cold 
front is detected also in the MOS image. More importantly 
another sharp feature (considering also the XMM PSF at this 
off-axis angle, ^ 50" 90% encircled energy fraction at 1.5 
keV) is clearly revealed at ^350", i.e. 65 kpc (see right panel 
of Fig. 6). The feature is detected also in the pn data. 

4. SPECTRAL ANALYSIS AND TEMPERATURE MAPS 

Details about the Chandra and XMM spectral extraction 
and analysis can be found in Humphrey et al. (2006) and 
Gastaldello et al. (2007b). Here we briefly summarize that for 
Chandra appropriate count-weighted spectral response matri- 
ces were generated for each region using the standard CIAO 
tasks mkwarf and mkacisrmf and for XMM using the 
SAS tasks rmfgen and arf gen in extended source mode. 



For each spectrum, we estimated the background through lo- 
cal modeling using the method outlined in Humphrey et al. 
(2006) and Gastaldello et al. (2007b). The spectra were re- 
binned to ensure a S/N of at least 3 and a minimum 20 counts 
bin"'. 

4. 1 . The inner 10 kpc 

We extracted spectra from the Chandra dataset in a series 
of regions as depicted in Fig. 7, determined using the con- 
tour binning technique (Sanders 2006) with a S/N=50. We 
fitted the background subtracted spectra with an APEC ther- 
mal plasma model (Smith et al. 2001) with the absorbing col- 
umn density fixed at the Galactic value (Dickey & Lockman 
1990). Solar abundances are in the units of Grevesse & Sauval 
(1998). To account for the undetected point-sources we added 
a 7.3 keV bremsstrahlung component (this model gives a good 
fit to the spectrum of the detected sources in nearby galax- 
ies, Irwin et al. 2003) for all the regions, because they fall 
within the twenty-fifth magnitude isophote (D25) of the cen- 
tral galaxy NGC 5044. 

The temperature map thus obtained is shown in Fig. 7. The 
clearest feature is the cool gas present in the SW (Sector III) 
arm (cooler than any other region in the map) which is spa- 
tially coincident with the Ha and dust filament (Temi et al. 
2007b, see Fig. 11). The SE arm (Sector I) is hotter. There is 
presence of cooler gas with respect to the surroundings in the 
direction of the inner surface brightness edge (i.e. region 26 
in Fig. 7). 

One possible concern about the interpretation of the spec- 
tral fits and consequently of the surface brightness features 
is the suggestion of limited multiphase gas in the tempera- 
ture structure of NGC 5044, as derived by the preference of 
two temperature (2T) models over single temperature (IT) 
models in radial annuli (Buote et al. 2003a). This is further 
complicated by projection effects, given the fact that we are 
dealing only with projected spectra in the analysis presented 
here: a proper deprojection is problematic due to uncertainty 
in the projection geometry and it is beyond the scope of the 
present paper. However it is evident in the radial analysis of 
Buote et al. (2003a) that the cool component clearly domi- 
nates in the inner 10 kpc (see their Figure 6) and that a single 
phase description, modulo projection effects, is an appropri- 
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Fig. 6. — Left panel: Surface brightness profiles in selected angular sectors, with PA measured from the N direction. The inner edge can be easily seen, though 
not as sharp as in the Chandra profile because of the XMM larger PSF. 

Right panel: Zoom over the interesting radial region for the outer edge. The characteristic shape of a rapid change in the slope of the profile across the edge can 
easily be seen in the sector PA 280-350. The profile in the undisturbed sector PA 350-80 is plotted for comparison. 



ate description of the data for these inner regions. We inves- 
tigated fitting 2T models as in Buote et al. (2003a) with the 
addition of the bremsstrahlung component for the unresolved 
point sources. The results of the two-dimensional analysis are 
consistent with the radial analysis of (Buote et al. 2003a) and 
with the dominance in terms of emission measure of the cool 
component over the hot one. 

To investigate in more detail the nature of the bright soft 
X-ray filaments we extracted spectra representative of the two 
southern filaments in two 30"x45" rectangular regions at a 
distance of 30" from the center (the blue boxes shown in 
Fig. 7); the spectrum taken from the SW filament has a Fe- 
L feature with lower excitation energy lines more prominent 
than the spectrum from the SB filament, signature of a lower 
temperature (see Fig. 8). If we fit the two spectra with a IT 
model we obtain a temperature difference which is significant 
at 9cr (kT = 0.65 ±0.01 keV for the SB arm and 0.78 ±0.01 
keV for the SW arm). The regions connected to the filaments 
are cooler compared to all the other regions in the map (by 
3.6(7 compared to the central region #1, which has the second 
lowest temperature in the map, 0.70 ± 0.01). If we try to fit 
with a 2T model the spectrum of the SB arm we obtain only a 
small decrease in ~x^ (x^/dof = 57/40 for the 2T model com- 
pared to 60/42 of the IT model). If we repeat the exercise for 
the regions of brighter emission surrounding the northern cav- 
ity (the cyan boxes shown in Fig. 7 encompassing the northern 
filaments highlighted in Fig. 1) we obtain 0.74 ± 0.02 keV for 
the NW filament and 0.82 ± 0.02 keV for the NB filament. 

4.2. The surface brightness edges 

To investigate the nature of the surface brightness edges we 
have to determine the temperature across the surface bright- 
ness jump. We therefore extracted Chandra spectra for the 
inner jump together with "control' spectra taken from regions 
at the same radial distance from the center but at different po- 
sition angles. We applied the same method to the outer edge 
when extracting XMM spectra. In Fig. 9 we plot the tempera- 
ture across the edges obtained by fitting a IT model. In Table 
2 we present the results of IT and 2T fits for selected regions 
inside and outside the edges together with control regions at 
the same radial range but at different position angles. 



For the inner edge a clear temperature jump is detected, 
with the brighter regions inside the edge cooler than the outer 
part. 2T models are preferred over IT models which points to 
some degree of projection effects in the temperature determi- 
nation across the edge, probably caused by the cool gas which 
is producing the edge. Other regions covering the same range 
of radii but at different position angles do not show in-fact the 
same behavior, with the exception to some degree of the re- 
gion immediately close to the edge. We extracted XMM MOS 
and pn spectra from the same regions inside and outside the 
inner edge finding consistent results within the la errors for 
both IT and 2T models. 

For the outer edge the temperature jump is more subtle and 
of low statistical significance (1.6cr), but, contrasted with the 
slightly declining trend of the control regions (we are in the 
region around the peak of the azimuthally averaged temper- 
ature profile), it is consistent with the picture of cooler gas 
inside the edge. 

5. SYSTBMATIC BRRORS 

This section contains an investigation of possible system- 
atic errors in the data analysis relevant for this paper A thor- 
ough analysis of many of the issues involved has been con- 
ducted in Buote et al. (2003a) and Buote et al. (2003b). 

5.1. Image Binning 

The presence of the morphological structures discussed in 
the paper, in particular the disturbed morphology with fila- 
mentary structure of the inner regions, is already evident from 
the raw Chandra image. The use of the particular binning 
technique of Sanders (2006) does not introduce any spuri- 
ous feature. We also used the Weighted Voronoi Tessellation 
(WVT) binning algorithm of Diehl & Statler (2006), which is 
more robust against the introduction of spurious features than 
adaptive smoothing, as provided for example by the CIAO 
task csmooth and we found the same features. 

5.2. Plasma Codes 

We compared the results obtained using the APBC code 
to those obtained using the MBKAL code (Kaastra & Mewe 
1993; Liedahl et al. 1995) to assess the importance of different 
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Fig. 7. — Left panel: Temperature map of the core regions of NGC 5044 based on the SN=50 contour binning. The map has been smoothed for better 
presentation with a Gaussian kernel of 6". Color bar units are in keV. The magenta cross marks the center of NGC 5044. The spectra presented in Fig. 8 and 
discussed in the text have been extracted from the rectangular regions shown in blue and they encompass the SE and SW filament; the spectra extracted from the 
cyan boxes and discussed in the text encompass the NE and NW filaments highlighted in Fig.l;. Right panel: Values and Icr error bars for the temperature in 
each region of the temperature map. 
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energy (keV) 

Fig. 8. — The spectra from the rectangular boxes taken respectively from 
the SW arm (black data points) and the SE arm (red data points). Note the 
shape of the Fe-L "hump" and the prominence of lower energy lines, signa- 
ture of lower temperature, in the spectrum of the SW arm. 

implementations of the atomic physics and different emission 
line lists in the plasma codes. We found no qualitative dif- 
ferences between the two codes: the fitted temperatures agree 
to within 5% and abundances within 10-20% (where IT and 
2T models are compared accordingly) and well within their 
statistical errors. 

5.3. Bandwidth 

We explored the sensitivity of our results to our default 
lower limit of the bandpass, /s^in = 0.5 keV. For comparison 
we performed spectral fits with £■„!„ = 0.4 keV and E^^n = 
0.7 keV. The fitted temperatures are consistent between mod- 
els; iron abundances are better constrained, in particular for 
2T models, when using the larger bandwidth (see Buote et al. 



2003b). 



5.4. Variable Nu 



We take into account possible deviations for A^h from the 
value of Dickey & Lockman (1990) allowing the parameter 
to vary by ±25%. We found no qualitative differences be- 
tween the two cases. If we leave A^h free to fit we derive 
higher values in the core, but still consistent at Ict with the 
Galactic value in the majority of the regions of Fig. 7. All 
the other spectral parameters are basically unchanged. This 
is not due to excess absorption but just the exploitation of an 
additional free parameter by the fitting program to model de- 
viations from a single temperature model. 

5.5. Background 

Since NGC 5044 is sufficiently bright and the temperature 
determination stems from the measurement of the shape and 
width of the Fe-L shell which is a prominent feature in the 
spectra, the fitted temperature values are quite insensitive to 
errors in the background normalization. 

6. DISCUSSION 

6.1. The nature of the cavities and filaments in the inner 10 

kpc 

The X-ray analysis shown in this paper reinforces the sce- 
nario proposed for NGC 5044 by Temi et al. (2007b), which 
we briefly summarize below. In a recent survey of ellipti- 
cal galaxies observed with the Spitzer telescope Temi et al. 
(2007a) found spatially extended cold interstellar dust emit- 
ting at 70 /im around many group-centered X-ray luminous 
elliptical galaxies. The source of this dust is the dust-rich 
cores of the central galaxy, not from mergers (Mathews & 
Brighenti 2003; Temi et al. 2007a). Since this dust is thought 
to be in direct contact with the hot, virialized interstellar gas 
(kT ~ 1 keV), it has a short lifetime (~ 10^ yrs) to sputter- 
ing destruction by thermal ions. Consequently, this dust is a 
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Fig. 9. — Left panel: Values and \cr error bars for the temperature across the inner surface brightness edge and in control sectors. 
Right panel: Values and \a en'or bars for the temperature across the outer suiface brightness edge and in control sectors. 
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Parameters from the spectral fits across the edges 
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Note. — Results of IT and 2T spectral fits for selected regions as described in the text in §4.2. The first and second columns refer to the radial range (in 
arcseconds) and azimuthal range (PA in degrees) of the fitted sectors discussed in §4.2. The norm parameter is the emission measure of the APEC model as 
defined in Xspec: I0^^'*(f neni,dV)/4TTD^(l +z)^ with units cm"^. Tcand Threfer to the temperature of the cool and hot component of the 2T model, whereas 
normc and norm;, refer to the con'espondent emission measures. No entry for 2T models means that there was no improvement over a IT model. Regions with 
no entry for Ti, did not significantly constrain that parameter 



spatial tracer of extremely transient events that recently oc- 
curred on kpc scales. NGC5044, together with NGC 4636, 
has been investigated in more detail in Temi et al. (2007b). 
NGC 5044 is not only extended at 70 /im but also at 8 /im in a 
manner similar to the highly asymmetric Ha optical line emis- 
sion from warm gas. Current evidence is consistent with the 
hypothesis that dust has been buoyantly transported from the 
galactic cores out to several kpc following a feedback heating 
event. Furthermore, disorganized fragments of optically ab- 
sorbing dusty gas are visible in HST images within the central 
'^100 parsecs, providing further evidence of a recent central 
energy release. In particular, the association of interstellar 
PAH emission (responsible of the 8 /im emission) and warm 
gas {T ~ lO'* K) in NGC 5044 indicates that we may be view- 
ing this galaxy at a rare moment immediately following a re- 
lease of energy near the central black hole. 



The Chandra data show the impact of this energy release 
on the hot X-ray emitting gas which is undergoing short-term 
gas-dynamical activity. The detection of a pair of X-ray cav- 
ities allow us to give an estimate of the amount of energy 
deposited from the AGN and the time of its release, follow- 
ing a now standard analysis (e.g. Birzan et al. 2004; Dunn & 
Fabian 2006). For each of the two cavities a size and posi- 
tion was measured, assuming the projected shape is an ellipse 
measured from the exposure-corrected, un-processed image 
as depicted in Fig. 10. We take as a guide the cavity on the 
south because of the higher contrast and the surface brightness 
profiles of Fig. 3. The two cavities are at a projected distance 
/? of 33" (6.1 kpc) from the center of the bubble to the center 
of NGC 5044, with a projected semi-major axis a of 17" (3.1 
kpc) and a projected semi-minor axis /? of 9" ( 1 .7 kpc). We do 
not consider here the dip in surface brightness at PA 212-244 
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TABLE 3 
Cavities time scales and derived powers. 
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0.41 
0.72 


32.1 ±0.09 



"It refers to the sum of the pV energies of the two cavities 
The first value refer to g calculated according to the hydrostatic equilibrium mass estimate whereas the second entry refer to the 
estimate made using the stellar velocity dispersion of the central galaxy, as described in the text. 

'^The power is the pV /t work only, with 7/(7- I) not accounted for For a fully relativistic plasma the values for the powers 
need to be multiplied by 4. 

Bolometric (0.1-100 keV) X-ray luminosity within a radius of 27 kpc, where the cooling time is 3 Gyr. 




Fig. 10. — Exposure corrected 0.5-5 keV image of NGC5044. Elliptical 
regions in black represent the shape and size of the detected X-ray cavities. 



because it not a clear deviation from the azimuthally averaged 
surface brightness profile. 

Assuming that the cavities are bubbles completely devoid of 
gas at the local ambient temperature and taking the pressure 
and temperature of the gas surrounding the cavities to be the 
azimuthally averaged values at the projected radius of its cen- 
ter as calculated in the analysis of Gastaldello et al. (2007b), 
we can calculate the PdV work done on the X-ray gas by the 
AGN in inflating the cavities. The actual third axis c of the 
ellipsoid is unknown: we assume a prolate shape (b = c). In 
this case we obtain pV = 2.29 ± 0.08 x 10^^ ergs for each cav- 
ity (the gas pressure at 6.1 kpc is 2.08 ±0.07 x 10~" dyne 
cm"^). Following Birzan et al. (2004) we define three dif- 
ferent time scales and therefore potential ages for the bub- 
bles: the sonic time scale ff^ = R/c,, where c, = 445 ± 3 km/s 
taking kT = 0.77 ±0.01 at 6.1 kpc; the buoyancy time scale 

thuoy =R/ ySC/lgV where V is the volume of the bubble, S is 
the cross section of the bubble and C = 0.75 is the drag coef- 
ficient (Churazov et al. 2001). The gravitational acceleration 
g can be calculated either as ,g = GM(< R)/R^ using the mass 
profile derived in Gastaldello et al. (2007b) or, as there can 
be reasonable doubts about the strict applicability of hydro- 
static equilibrium in the inner kpc of NGC 5044 (we are likely 
underestimating mass and therefore overestimating tbuoy, as 
we argued in Gastaldello et al. 2007b, given also the failure 
to detect the stellar mass of the central galaxy), we use the 



stellar velocity dispersion of the central galaxy, under the ap- 
proximation of an isothermal sphere and calculate g « 2(t^ /R 
(Binney & Tremaine 1987); the time required to refill the dis- 
placed volume as the bubble rises upward, tr = 2y/ r/g, again 
using for g either the X-ray determined mass or the stellar 
dispersion velocity. The calculated order-of-magnitude esti- 
mates for the time scales described above are shown in Table 
3. and they are consistent with the short sputtering lifetime for 
the extended dust, ^10^ yrs, discussed in Temi et al. (2007b) 
and point towards the same episode of AGN feedback. 

We can compare estimates of the energy released by the 
AGN to the radiative losses within the cooling region, here 
estimated to be the region within a radius where the cooling 
time is less than 3 Gyr (e.g. Dunn & Fabian 2006), also tab- 
ulated in Table 3 which corresponds to 27 kpc. The task of 
deriving the outburst energy from the observable pV is com- 
plicated by various unknowns, like the adiabatic index 7 of 
the material inside the bubble (7 = 5/3 if thermal or 7 = 4/3 
if non-thermal, see for example the discussion in Mathews 
& Brighenti 2008) or the correct hydrodynamic evolution of 
the bubble expansion, like its initial over-pressure compared 
to the ICM or the possibility of continuous inflation, to ex- 
plain the paradoxical result of bubble energy increasing with 
distance from the center of the gravitational potential (Diehl 
et al. 2008). Using the commonly adopted recipe of estimat- 
ing as 4pV, with 7 = 4/3, the energy input arising from the 
cavities alone, this energy can be at most ~ 8 x 10"*' erg s"', 
which falls short to counteract cooling losses. 

It is remarkable that a pair of cavities close to the nuclear 
source are lacking extended high frequency radio emission. 
This is contrary to what is generally observed in particular in 
clusters of galaxies. Analysis of VLA observations at 1.56 
and 4.9 GHz confirms the presence of just a point source 
with beam-deconvolved size < 3"and with a very flat (a ^ 0) 
spectral index (Giacintucci et al., in preparation). Two other 
sources sharing the same behavior are NGC 4636 (e.g., Jones 
et al. 2002) and NGC 4552 (Machacek et al. 2006). However 
the cavities in these sources have been interpreted as caused 
by a nuclear outburst which is directly causing shock heat- 
ing, also because the brighter rims surrounding the cavities 
are actually hotter than the surrounding medium, rather than 
displacement by the radio lobes (Jones et al. 2002). Or the 
cavities could just be related to a previous outburst and might 
just be ghost cavities, as suggested by O' Sullivan et al. (2005) 
for NGC 4636 and the same reasoning could apply to NGC 
5044. In that case it remains to be explained why the syn- 
chrotron emission has faded to undetectable values in these 
sources and not in cavities with approximately the same ages 
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Fig. 1 1 . — Left Panel: Ha + [N ii ] map of NGC5044 taken from Caon et al. (2000). Right Panel: Surface brightness contours of the Ha + [N ii] map 
superimposed on the contour binned image of Fig. 1 . 



in many other sources (e.g., Birzan et al. 2004). Other ghost 
cavity systems with a weak central radio source bearing no 
obvious relation to the observable cavities which are close to 
the central galaxy are HCG 62 (Morita et al. 2006) and NGC 
741 (Jetha et al. 2008, where only one cavity is present). Sen- 
sitive low-frequency radio observations will shed further light 
on this issue. 

The association of dust. Ha and soft X-ray emission 
showed by Temi et al. (2007b) is also strengthened by the 
analysis in this paper The presence of the N cavity and in 
particular of the relatively cooler emission in the NW filament 
(see §4.1) likely explains the origin of the N Ha filament. X- 
ray filaments are present at both sides of the cavities, but only 
the ones with the presence of dust are showing optical emis- 
sion and cooler X-ray emission (see Fig. 11 and the analysis 
of §4.1). Temi et al. (2007b) showed that the cospatiality of 
these features can be explained as the result of dust-assisted 
cooling in an outflowing plume of hot dusty gas: dust can cool 
buoyant gas to 10"^ K, which emits the optical emission lines 
observed. The warm-gas phase is maintained in thermal equi- 
librium near ~ 10"* K by radiative losses and likely a combi- 
nation of thermal conduction and UV heating from post-AGB 
stars (see the calculation of Tab. 5 in Macchetto et al. 1996, 
ionization from post-AGB stars could explain only 42% of 
the optical line emission in NGC 5044). It is unfortunate that 
the deeper Ha observation of Caon et al. (2000), compared 
to the one by Goudfrooij et al. (1994b) used in the compar- 
ison of Temi et al. (2007b), is affected by a CCD defect in 
the southern region co-spatial to the X-ray cavity: it looks 
like the nebular emission is more extended and can cover also 
the cavity. Another tantalizing evidence of the association of 
these three components comes for example from NGC 5846 
(Trinchieri & Goudfrooij 2002, and references therein). Evi- 
dence of PAH emission in Spitzer IRS spectra has been found 
also in NGC 1275 (Johnstone et al. 2007). It has been dis- 
cussed that AGN feedback can provide some heating (nega- 
tive feedback), but it can also be responsible for a positive 
feedback, i.e. for inducing cold gas production (e.g., Pizzo- 
lato & Soker 2005; Revaz et al. 2008). In this scenario AGN 
feedback itself is responsible for the production of overdense 
blobs which cool rapidly producing Ha emitting gas, molec- 



ular gas and star formation. Dust transported from the dusty 
disks of the central elliptical into the ICM by the episode of 
AGN feedback acts as a catalyst for the cooling of the gas. 
It would be therefore interesting for the recent renewed the- 
oretical interest in the generation and survival of the optical 
line emission filaments (Nipoti & Binney 2004; Revaz et al. 
2008; Pope et al. 2008) to include the neglected dust-assisted 
cooling in the energy balance equation. 

6.2. The cold fronts and the dynamical state of the group 
NGC 5044 

The spectral analysis presented in §4.2 indicated that the 
pair of surface brightness discontinuities detected in the 
Chandra and XMM images are cold fronts. Given the large 
scale relaxed morphology of NGC 5044 and its rising entropy 
profile (Gastaldello et al. 2007a) as commonly observed for 
bright relaxed groups, the source nicely fits in the explanation 
for the emergence of these features proposed by Ascasibar & 
Markevitch (2006) as due to gas sloshing caused by an off- 
axis merger with a smaller satellite. Another tell-tale sign 
of this encounter, as also suggested by Ascasibar & Marke- 
vitch (2006), is the presence of a peculiar velocity of the 
central galaxy, which is often exhibited by central cD galax- 
ies in relaxed clusters (e.g., Oegerle & Hill 2001). Indeed 
NGC 5044 is known to have a ~ 150 km/s peculiar velocity 
with respect to the mean group velocity (Mendel et al. 2008, 
and references therein). A sub-clump of galaxies statistically 
significant to the Dressler-Schectman test (Dressier & Sheet- 
man 1988) has been detected in the north-east outskirts of the 
NGC 5044 (Mendel et al. 2008) which might be the smaller 
merging sub-group responsible for the gas sloshing and the 
cold fronts. For a virial mass for NGC 5044 of 3.7 x 10'^ M© 
(Gastaldello et al. 2007b), assuming a mass ratio of 5 as in 
the reference case of Ascasibar & Markevitch (2006) would 
imply a mass of ~ 8 x 10^^ M© for the satellite. If this can- 
didate sub-clump is real, it would corroborate the proposed 
scenario and the mechanism proposed by Ascasibar & Marke- 
vitch (2006) for the formation of cold fronts (see also the case 
of the cluster A 496, Dupke et al. 2007) 

Surface brightness discontinuities as cold fronts have so far 
been discussed in many merging elliptical galaxies and groups 
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(e.g., NGC 1404, Machacek et al. 2005). An unusual discon- 
tinuity for the group NGC 507, interpreted as an abundance 
jump, closely related to the expansion of a radio lobe, has 
been discussed by Kraft et al. (2004). NGC 5044 is to our 
knowledge the first relaxed group for which cold fronts have 
been discussed in close similarity to the ubiquitous ones de- 
tected in relaxed clusters. Further examples, given enough 
data quality, are likely to be discovered: suggestions have 
already been made for objects like MKW 4 and IC 1860 
(Gastaldello et al. 2007b). 



The detection of such a rich phenomenology has been pos- 
sible due to NGC 5044 being bright and nearby, much like on 
bigger mass scales has been possible with M87 and Perseus. 
As for these objects up-coming deeper X-ray observations 
with Chandra (David 2007) and XMM (Kaastra 2007) will 
shed further light on the properties of NGC 5044 and groups 
of galaxies in general. Multi-wavelength observations, i.e. in 
Ha and radio bands, are needed for a deeper understanding of 
the physical processes in the core. 



7. CONCLUSIONS 

We have presented results of a two-dimensional analysis of 
the currently available Chandra and XMM data for the bright 
nearby galaxy group NGC 5044. The results can be summa- 
rized as follows: 

- A pair of X-ray cavities have been detected, further con- 
firming the recent outburst indicated by the extent and 
morphology of Ha and dust emission. 

- The presence of cooler filamentary X-ray emission co- 
spatial with Ha and dust emission strengthens our pre- 
vious suggestion that dust-aided cooling could con- 
tribute to the Ha emission. 

- The presence of a set of two cold fronts together with a 
peculiar velocity of the central galaxy NGC 5044 sug- 
gests a disturbance of an overall relaxed system by an 
off-axis merger with a smaller satellite. 
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